Existence of the excitonic insulator phase in the extended Fahcov-Kimball model: an 

S'0(2)-invariant slave-boson approach 
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We re-examine the three-dimensional spinless Falicov-Kimball model with dispersive / electrons 
at half-filling, addressing the dispute about the formation of an excitonic condensate, which is 
closely related to the problem of electronic ferroelectricity. To this end, we work out a slave-boson 
functional integral representation of the suchlike extended Falicov-Kimball model that preserves 
the S0{2) (g) U{1)®'^ invariance of the action. We find a spontaneous pairing of c electrons with / 
holes, building an excitonic insulator state at low temperatures, also for the case of initially non- 
degenerate orbitals. This is in contrast to recent predictions of scalar slave-boson mean-field theory 
but corroborates previous Hartree-Fock and RPA results. Our more precise treatment of correlation 
effects, however, leads to a substantial reduction of the critical temperature. The different behavior 
of the partial densities of states in the weak and strong inter-orbital Coulomb interaction regimes 
supports a BCS-BEC transition scenario. 

PACS numbers: 71.28.-(-d, 71.35.-y, 71.35.Lk, 71.30.-fh, 71.2S.+A. 71.27. +& 
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I. INTRODUCTION 

The excitonic instability in solids is driven by the 
Coulomb attraction between electrons and holes which 
under certain conditions causes them to form bound 
states. At the semimetal-semiconductor transition the 
conventional ground state of the crystal may become un- 
stable with respect to the spontaneous formation of exci- 
tons. Starting from a semimetal with a sufficiently small 
number of electrons and holes, such that the Coulomb 
interaction is basically unscreened, the number of free 
carriers will vary discontinuously under an applied per- 
turbation^, signaling a phase transition. Approaching 
the transition from the semiconductor side, an anomaly 
occurs when the (indirect) band gap, tuned, e.g., by ex- 
ternal pressure, becomes less then the exciton binding 
energy^. As a consequence, a new distorted phase of the 
crystal, with spontaneous coherence between conduction 
and valence bands and a gap for charged excitations, de- 
velops. It separates, below a critical temperature, the 
semimetal from the semiconductor. This state, known as 
'excitonic insulator' (EI), can be regarded as an electron- 
hole pair (exciton) condensate^. By nature, depending on 
from which side of the semimetal-semiconductor transi- 
tion the EI is approached, the EI typifies either as a BCS 
condensate of loosely bound electron-hole pairs or as a 
Bose-Einstein condensate (BEC) of preformed tightly- 
bound excitonsi"— . 

The mean-field description of the EI phase is very sim- 
ilar to the BCS theory of superconductivity, and has 
been worked out long-time a.gcS^—. In this context 
also the transition from BCS to BE condensation was 
discussed^ii^^— . Surprisingly enough, the quantitative 
semimetal-EI-semiconductor phase diagram has been de- 
termined only quite recently^"—. All these investiga- 
tions, having normal (excited) semiconductor systems in 
mind, rest on the standard effective-mass Mott-Wannier- 



type exciton model. Thereby important band-structure 
and correlation effects, as well as the exciton-excition and 
exciton-phonon interaction, and the inter-valley scatter- 
ing of excitons were largely neglected. 

In nature, the EI state is evidently rare. One obsta- 
cle for creating an excitonic condensate is the far-off- 
equilibrium situation caused by optical excitation of ex- 
citons. But also in thermal equilibrium situations, EI 
states are expected to occur only under very particu- 
lar circumstances, e.g., if conduction and valence bands 
are adequately nested^^. Actual materials experimen- 
tally studied from the viewpoint of the EI are numbered. 
One recent example is quasi-one-dimensional Ta2NiSe5 
with highly polarizable Se, where angle-resolved photoe- 
mission spectra reveal an extreme valence-band top flat- 
tening indicating that the ground state might be viewed 
as an EI^^. At present, the transition-metal dichalco- 
genide ir-TiSe2 seems to be the only candidate for a 
low-temperature phase transition to the EI without the 
influence of any external parameters other than the tem- 
perature. Here the onset of an EI phase was invoked as 
driving force for the charge-density-wave transitioi*22i2i. 
Semiconducting, pressure-sensitive mixed valence mate- 
rials, such as TmSeo.45Teo.55, are further candidates for 
exciton condensation. Fine tuning the excitonic level to 
the narrow 4/ valence band, a rather large number of 
about 10^*^ — 10^^ cm~'^ (small-to-intermediate sized) ex- 
citons can be created, which presumably condense into 
an EI state, at temperatures below 20 K in the pres- 
sure range between 5 and 11 kbar^^"^^. Clearly, for these 
rather complex transition- metal/rare-earth compounds 
with strong electronic correlations, simple effective-mass- 
model based theories will be too crude. 

The investigation of Falicov-Kimball-type models of- 
fers another promising route towards the theoretical 
description of the EI scenario. In its original form, 
the Falicov-Kimball mode P^'^^ contains two types of 
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fermions: localized / electrons and itinerant c electrons 
with orbital energies Ef and , respectively. An on- 
site Coulomb interaction U of both species determines 
the distribution of the electrons between the / and c 
sub-systems, and therefore may drive a valence transi- 
tion, provided there is a way to establish f-c coherence. 
At first glance this can only be achieved by including 
a hybridization of the / and c bands-^'^*. It has been 
shown, however, that a finite / bandwidth, being cer- 
tainly more realistic than entirely localized / electrons, 
will also induce f-c transitions22i2£. The model with di- 
rect c-c and /-/ hopping (oc ^c//), is sometimes called 
extended Falicov- Kimball model (EFKM), the Hamilto- 
nian of which takes the form 

H =E, + E + ^/ E + E 
+ c^E (1) 

i 

Here // (cj) creates an / (c) electron at lattice site i, and 

'^i/ — fifi — c|c,j) are the corresponding number 
operators. Let us emphasize that the / and c bands 
involved have different parity^S. For tftc < (t/ic > 0), 
we may have a direct (indirect) band gap. For tf = 
(dispersionless / band), the local / electron number is 
strictly conserved^. 

In the past few years, both the Falicov-Kimball model 
with hybridizatioi*2L28^-^ and the EFKM^s-^iSS ^^^^ 

been studied in connection with the exciting idea of 
electronic ferroelectricity. The origin of electronic fer- 
roelectricity is a spontaneously broken symmetry due 
to a non- vanishing (c^/) average, which causes finite 
electrical polarizability without an external, interband- 
transition driving field. As (c^/) is basically an ex- 
citonic expectation value, indicating the pairing of c 
electrons with / holes, the problem of electronic ferro- 
electricity is intimately connected with the appearance 
of an excitonic condensate. Accordingly, the question 
whether the ground-state phase diagram of the EFKM 
exhibits an EI state has attracted much attention. By 
means of constrained path Monte Carlo (CPMC) tech- 
niques the T — phase diagram of the EFKM was de- 
termined in one and two dimensions in the strong and 
intermediate-coupling regime o^^i'^° . In both cases a ferro- 
electric phase was detected. A subsequent Hartree-Fock 
calculation shows that the mean-field phase diagram of 
the two-dimensional EFKM agrees even quantitatively 
with the CPMC data^-, supporting the applicability of 
Hartree-Fock and RPA schemes to three- and infinite- 
dimensional systems^i^i^. Surprisingly, the more so- 
phisticated scalar slave-boson theory failed to find the EI 
phase when the / and c orbitals are non-degenerate^I. 

The continued controversy, regarding the existence of 
the EI phase in the EFKM, motivates us to re-examine 
the problem using an improved auxiliary boson approach 
that ensures the rotational and gauge symmetries of the 
EFKM within a functional integral scheme. 



II. SLAVE-BOSON THEORY 

A. Slave-boson Hamiltonian 

The extended Falicov-Kimball Hamiltonian ([T|) can be 
rewritten as an asymmetric Hubbard model if the or- 
bital flavor (/, c) is represented by a pseudo-spin variable 
(cr =t,i)- Using the spinor representation 

where the vectors (a|^) are built up by the fermion 
annihilation (creation) operators a^j^ = fj^^^ and af^ = 
cl'\ H becomes 

H^J2 + E ^ + ^ E '^^'^^T • (3) 

In Eq. (12), K — t^/ti gives the ratio of the /- and c- 
bandwidths (t^ = 1 fixes the unit of energy). Obviously, 
the usual Hubbard model follows for E^ — E^ and k = 1 . 
Without loss of generality we choose E^ — Q and E^ < 
in what follows. 

Now the slave-boson representation of the EFKM is 
constructed by replacing the fermionic Hilbert space by 
an enlarged one of pseudo- fermionic and bosonic states. 
The local states, representing the original physical states 
of the EFKM ([T|) in the enlarged Hilbert space in a one- 
to-one manner, can be created in the following way: 

|0.) ^ e^lvac), (4) 
1 2,) &l^&yi\Yac), (5) 

W^) ^ E^l^'Ul^^^)- (6) 
p 

The pseudo-fermions dip satisfy anti-commutation rules 
{fljp, a]p/} = SijSppi, while usual Bose commutation rules 
hold for the slave bosons: [e,,e]] ^ [p^p,p^,p]p,pj = 
5<5y <^piP4<5p2p3 , and [d^jd]] ~ 6ij- The boson number op- 
erators, ejcj, 2Trp|p., and d\d^ project on an empty, 
a singly occupied and a doubly occupied state, respec- 
tively. By introducing a slave-boson matrix-operator for 
the case of single occupancy, p'^^\ we adapt the spin- 
rotation-invariant slave-boson formulation of the Hub- 
bard model2S (for a generalization to multi-orbital mod- 
els see2S), in order to avoid difficulties that may arise 
from the scalar nature of the bosons in approximative 
treatments22ii£. The decomposition 

p=0 

into scalar (singlet) p|J'' and vector (triplet) p^^^ = 
iP^ix ^Piv ^P^iz) components, where Tq is the unity matrix 
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and r the vector of Pauli spin matrices, is given as 



(t) _l / P^o +P 



(t)^Jt) (t) 



(t)> 



VP 



(t) 



ip. 



(t) „(t) 



P 



(t) 



(8) 



Of course, it is crucial to select out of the extended 
fermion-boson Fock space the physical states. This can 
be achieved by imposing two sets of local constraints: 



(1) 



(2) 



ele, +2Trptp^ -1 =0, (9) 

a, ® aj + 2 ptp . + dld^ To - To = . (10) 



Cp'' expresses the completeness of the bosonic operators, 
i.e. each lattice site i can only be occupied by exactly 

(2) 

one boson. C_l relates the pseudo-fermion number to 
the number of p and d bosons. 

Correspondingly, the mapping of the physical electron 
operators into products of new pseudo-fermions and slave 
bosons in the hopping term of H is 



3.4 } a,; . 



The choice of the hopping operators Zj is not unique. 
This arbitrariness can be used, e.g., to reproduce, for 
the Hubbard model case, the correct free-fermion result 
at [/ = and the Gutzwiller result for any finite V at 
the mean-field level, where the constraints ^ and 
are fulfilled only on the average. This is guaranteed by 
choosing^Sj^i 

= L,e\M.,p^N, + L^M^d^N, (12) 

with 

= [(l-4d,)To-2ptpJ-i/2^ (13) 
- [(1 - 4e,)To - 2ptpJ-i/2 ^ (14) 
A-k - [l + ele, + 4d, + 2TrBt„.]i/2, (15) 



of H can be bosonized via 



and pip'p, = PP'Pi-p'-p ■ The Hubbard interaction term 



(16) 



That is, the transformation to slave-boson fields results 
in a linearization of the interaction, and we end up with 
the EFKM Hamiltonian in the form 

H = ^ a| (To + T, )a, -I- aki +uJ24d^■ 



B. Functional integral representation 



(17) 



To proceed further, it is convenient to represent the 
grand canonical partition function of the constrained sys- 
tem (fT7|) . Z = Tr e~^*^^~^^=\ as an imaginary-time path 



integral over Grassmann fermionic and complex bosonic 
fieldE42 



Z = 



D[Rp, -a^] D[e*,e] D[p;,p^] D[d*,d] d[X^'^ d[A|,2)] 



-/ drLir) 



X e 
with the Lagrangian 

L{r) 



(18) 



E 



+2Trpf ((9, + aW)to-A: 



+d*id, + \\'^ +U-TT^'>)d, 

+ J2 Hidr - M)ro + A1'^ + ^(To + T J) 



(2)^ 



Sii z* t Zj Sij . (19) 



Here = l/T is the inverse temperature and the time- 

(1) (2) (2) 

independent Lagrange multipliers , — Ko Zo + 

-'(2) -t , 

T are introduced to enforce the constraints via the 
integral representation of the 5-function^2.^ 

c+2Tii/l3 

= ^ J dAW e-^^'''^" (20) 

c 

(the path of the A-integration is parallel to the imaginary 
axis and one finds A G R+ at the physical saddle point). 

Next, exploiting the gauge symmetry of the action, we 
perform local time-dependent phase transformations: 



Gi r Ci 6 , 

d, ^ d, e-''^- , 

—i —i 

a, ^ a, e-*(^'«^o-v..zJ 



(21) 
(22) 
(23) 
(24) 



Note that both the original as well as the transformed 
Bose fields are complex. By the transformation (p4)) the 
kinetic contribution generates extra terms violating the 
50(2) (g) U{1)'^^ invariance of the model. Transforming 
the Lagrange multipliers into real time-dependent Bose 
fields, 



A, 



(1) 



(25) 



(26) 



and, in addition, restricting the phase transformation to 
50(2) ® C/(l)«^2 symmetry by 



(27) 
(28) 
(29) 
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the gauge invariance of the action is satisfied. We now 
make use of the gauge freedom to remove three phases of 
the Bose fields in radial gauge, where the fields are given 
as modulus times a phase factor: 



Si 

di 



(30) 
(31) 
(32) 



As a consequence, three bosons, e.g., ei(T), pio{T), 
and pizir), can be taken as real- valued, i.e., their ki- 
netic terms, being proportional to the time derivates in 
Eq. (jl9p . drop out due to the periodic boundary condi- 
tions imposed on Bose fields (0i(/3) = 0i(O)). However, 
the other three bosons pix, Piy, and di remain complex. 



di^\di\e with ijj, = ^j, - 2xjo + 



(33) 



„ , _, PiQ+Ptz {\Pix\ - i\piy\)e 



(34) 



This has to be contrasted to the 5/7(2) ® U{1) invariant 
Hubbard model {t-J model), where only one Bose field 
stays complex (all Bose fields become real)^"— . 
Using the familiar Grassman integration formula. 



D[dp,d]e 



, Trln[-G" 



(35) 



the grand canonical partition function can be represented 
as a functional integral over Bose fields only. 



Z = j D[e]D[p^]D[pl,p^]D[pl,Py]D[p,]D[d\d] 

Z?[A(i']D[A[,'V[A(2)]e-^ (36) 
with the effective bosonic action 

* ^ 

+ (Ap'+C/-2AL'')|d.p 
+P*ix9tPix +PiydTPiv + d*drd^ 
-Trln{-G^4,,,,(r,r')}, (37) 
where the inverse Green propagator is given by 



(2), 



"2^(ll0 + 21z)pp' -^i "^PP' 
~ {^.i tZ-j) pp' ,TT'i^ ~ ^ij) ■ 



S,^ d{T - r') 
(38) 



The trace in Eq. (1371) extends over time, space, and spin 
variables. 

The Hermitian matrix can be brought into the form 



where 





I 








1 


y^2) 





Pi 



■ Piz 



Pix ^Piy 
~Pi - Piz 



are the eigenvectors of p . , p. with 



Pi = |pd = y\Pix? 



■pIz^ 



(39) 

(40) 
(41) 

(42) 
(43) 



and 



{\-\d,\^)-\{p^^^P^f 



72 



(PiO + Pi) + '^i {PiO - Pi) 



Zi2 



(1 - e^) - i (PrO-Pif 
(l-|d,n-i(p,o-p,)^ 



(44) 



V2 



[Pia ~Pi) + di {piQ + Pi) 



il-ej)-^{p,o+P^)' 



Then we get 



^il2^jl(^'il ^12) : 
^ilViliZil ~ ^12) ■ 



(45) 

(46) 
(47) 



We note that only for the half-filled band case (n^ + 
1, ei = \di\), we find that zn — — Zi, i.e.. 



(48) 



becomes diagonal, and the matrix elements of the origi- 
nal Hamiltonian are reproduced by the slave-boson trans- 
formed model. That means, Eq. (p6|) with Eqs. (|37|) to 
(|48p provide an exact representation of the partition func- 
tion for the EFKM at half-filling. By contrast, for the 
5[/(2)-invariant Hubbard Hamiltonian with n — 1 and 
— the slave-boson representation of Z holds ex- 
actly for all fillings. 

For the EFKM case with t^t^ < (direct gap in the 
paraphase for large \E^\), the EI order parameter A_l and 
the 'Hartree shift' Az are respectively given ^.g^LS^j^ 



N 



t 



(50) 
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Using the constraints pO)) these relations can be ex- 
pressed as functional averages: 



U 



U 



(51) 
(52) 



C. Saddle-point approximation 

The evaluation of Eq. ([551) is usually carried out by 
a loop expansion of the collective action ([37]). At the 
first level of approximation, the bosonic fields are re- 
placed by their time-averaged values, and one looks 
for an extremum of the bosonized action with re- 
spect to the Bose and Lagrange multiplier fields (j)ia = 

„ - J \(2) J(2)\. 



dS 



= 



s^s 



The physically relevant saddle point {(f>ia} is deter- 
mined to give the lowest free energy (per site). 



/ = n/N + ^ln, 



(54) 



where, at given mean electron density n = + n^, the 
chemical potential fj, is fixed by the requirement 



1 on 



(55) 



0, — S//3 denotes the grand canonical potential. Clearly, 
an unrestricted minimization of the free energy is impos- 
sible for an infinite system, even within the static ap- 
proximation. Focusing on the possible existence of the 
EI phase, we consider only uniform solutions hereafter: 
{4'ia} = {4>a}- Note that the inclusion of a charge- 
density- wave phase is straightforward, e.g., by adapting 
the two-sublattice slave-boson treatment worked out for 
the Peier Is- Hubbard modet^^"'*®. 

Examining a tight-binding direct-band-gap situation 
in three dimensions, we have 



with 



and 



[l-<P-\{po+pY][l-dP-\{po~pY] 



7fc 



-2[cos kx + cos ky ~\- cos fc^] . 



(56) 



(57) 



(58) 



The trace in Eq. ([57)) can be easily performed in 
the momentum-frequency-domain after diagonalizing the 



propagator in pseudo-spin space. Then the free energy 
functional takes the form 

= X^^^e'+pl+p' + S-l) 

kv 

where 

n,,^ = [exp{/3(i?,,^ - A)} + l]-i (60) 
holds with the quasiparticle energies {v — ±.) 

= kjEt + i^+^yJk] (61) 
+iy^\[E^ + 2Ai'' + {k- 1)z27-]2 + (Af )2 . 



Here we have introduced A^'' = ±\/ {X^f')^ + {Xl"')^ , 



v(2)n2 , m(2)^ 



(53) P-L JpI+pI, and ji = X}; 



(2) 



Requiring that / becomes stationary with respect to 
the variation of the (j)a we obtain the following set of 
saddle-point equations: 



2po \ 2d? pI-p^ 



Ph-P 



,(2) ( Z^ 1 

^ 2po \2d^ pI-P^ 
1 1 

PaPz =2]^Z^^"^fc"£.' 



z^I, 



kv 



1 1 , r 



1 1 



Po=2 + 2V(l-^')(l-4pgp'): 



(62) 

(63) 
(64) 

(65) 

(66) 



d^^^[z\2-pl-p^) + 2pl 



- 2po^z^{2-pl-p^) + z^p^^pi\ , (67) 



U + 2A(^)^ + 2A(^)^ = '^" + ^^ zH 



^ Po ' Po 2pld2 



(68) 



with 



^ =(« + 1)^ Il7fenj, + ('* - H ^™fc7fc"fc. ' 

(69) 



kiy kv 



+ 2Xi^^ + (k - l)z^7j 



^{E^ + 2X^^^ + {n-l)z^j^r + {2X^l^y 



, (70) 



2X 



(2) 



^J{E^ + 2X^^^ + {n-l)z^j^)^ + {2X^l^y 



(71) 
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The EI order parameter (1491) and Hartree shift ((50)) be- 
come 

= UpoP± , (72) 
A, = 2UpoPz ■ (73) 

D. Zero temperature: BI to EI transition 

At zero temperature, the EFKM exhibits a trivial band 
insulator (BI) phase of a completely filled / (empty c) 
band {n ^ 1, < 0), provided the Hartree gap is finite: 



AniT = 0) = \E^ + 2A(2)| - 6(|k| + I) > . 



(74) 



That is to say, in the BI phase we have cP = 0, = 1, 
= 0, and in no way f-c coherence can develop: p± — 



(2) 



0. Then 



-I for all fc, and nj^_ 



1, = 



A 

resuh from Eq. with Eq. The constraint pO)) 

gives, together with Eq. ([64| . p^ — Pz = l/v^, leading 
to = 1, as for a non-interacting system. At the same 
time, / = according to Eq. ([M)) . and the correlation 

('2) 

equation (|68l) reduces to [/ + 2A2 = 0, which gives A^ — 
for the BI. 

Looking for an instability of the BI towards an EI state, 

(<2\ 

we find from Eq. ^ that AY^ = - Aj^ near the critical 
Coulomb interaction Uc2- Multiplying Eq. (|55|) by J7, for 
Aj^ 7^ 0, we get the T = gap equation 



1 V 

k 



Ai 



(75) 



which agrees with the Hartree-Fock result As a conse- 
quence, our 5'0(2)-invariant slave-boson approach repro- 
duces the BI-EI phase boundary of the EFKM Hartree- 
Fock ground-state phase diagra m^^i'^^ . At least for the 
2D case it has been demonstrated that this phase bound- 
ary agrees almost perfectly with that obtained by the 
CPMC method^i^. This also applies to our 5*0(2)- 
invariant slave-boson approach, e.g., for the 2D EFKM 
with K = —0.3 and U = 2, we obtain the critical value 
\E^^c\ = 3.23(3.26) for the EI-BI transition, using the 
2D tight-binding (square) density of states, in compari- 
son with ^ 3 29 (cf. Fig. 3 in Ref.lM). 

E. Scalar slave-boson approach 

If one contrariwise adopts the scalar slave-boson the- 
ory^ by introducing only four auxiliary bosonic fields 
per site, Ci, Pia, and d^, where pl^p.^^ (P4P4) projects 
on a singly occupied /- (c-) electron site i, the p. -matrix 
becomes diagonal: 



P. =2 



PiO + Pi. 







PiO - Pi. 



(76) 



(77) 
(78) 



That means, the 'spin- flip' terms in Eq. (ITOl) . 



do not occur. As these terms, in view of Eqs. (|49)) 
and (|5T|) . are essential for the formation of an excitonic 
insulator, the scalar slave-boson approach fails to de- 
scribe the EI phase, at least for finite orbital-energy dif- 
ference E^ ^ 0. At the (uniform) saddle-point level of 
approaximation, within scalar slave-boson theory, we find 
the band-renormalization factor 



(79) 



with p^ — V2ptt (Pi = V^Pii), and the correlation equa- 
tion (155)) simplifies to 



U = z^ 



PtP-l 



1 



(80) 



III. NUMERICAL RESULTS 

In the numerical evaluation of the self-consistency 
loop (IM1) - ([7T|) we proceed as follows: at given model pa- 
rameters Ef, K, and fixed total particle density n ~ 
rif + nc = 1, we solve the finite-temperature saddle-point 
equations for the slave-boson and Lagrange parameter 
fields {(j)a} together with the equation for the renormal- 
ized chemical potential /2 using an iteration technique. 
Thereby /c-summations were transformed into energy in- 
tegrals, introducing the (tight-binding) density of states 
for the simple cubic lattice. Convergence is assumed to 
be achieved if all quantities are determined with a relative 
error less than 10~^. Our numerical scheme allows for the 
investigation of different metastable states correspond- 
ing to local minima of the variational free-energy func- 
tional. Of course, we will always obtain a homogeneous, 
translational invariant solution without spontaneous (po- 
larization) exciton formation. In this case the / and c 

(2) 

bands are simply shifted by 2X1 , leading to a gapped 
band structure at large enough U. Besides these simple 
(semi-) metallic and BI phases, the T = Hartree-Fock 
ground-state phase diagram of the half-filled EFKM ex- 
hibits two symmetry-broken states^i^: the anticipated 
EI and a charge-density-wave phase. At Ef = Ec (degen- 
erate orbitals), the charge-density- wave ground state is 
stable for all values of k. It becomes rapidly suppressed, 
however, for Ef ^ Ec (non-degenerate orbitals), in par- 
ticular, if the c and / bandwidths are comparable'^^. As 
we are interested in the (uniform) EI phase only, we have 
confined our slave-boson approach to spatially uniform 
saddle points. With respect to charge-density-wave for- 
mation this will be uncritical for the parameter values 
studied in the following. 
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FIG. 1: (Color online) Stability region of the EI phase in the 
3D half-filled EFKM (the arrow marks the critical coupling 
where the Hartree gap ((74} opens). The inset shows the order 
parameter at zero temperature. Red dotted curves give the 
Hartree-Fock results for comparison. 





Figure [T] gives the slave-boson phase boundary of the 
EI in the U-T plane, calculated for = —2.4 and k — 
—0.8. Most notably, we obtain a stable EI solution for 
the non-degenerate band case, which has to be contrasted 
with the result of the scalar slave-boson approach^!. 

Let us first discuss the T = data. Here the numeri- 
cal semimetal-EI and EI-BI transition points at small and 
large Coulomb interaction, Uci — 0.74 and Uc2 — 9.3, re- 
spectively, agree with the Hartree-Fock results (see the 
dotted curve). The latter was proved analytically in 
Sec. Ill Dl The inset gives the ?7-dependence of the EI 
order parameter at T = 0. For Ud < U < Uc2, ^± 
only slightly deviates from the corresponding Hartree- 
Fock curve. 

The variation of the other bosonic fields is depicted in 
Fig. [2j where the solid curves belong to the parameter 
values used in Fig. [T] We see that the number of empty 
and double-occupied sites, and (P, is equal and goes to 
zero at the EI-BI transition, where we havep§ = pi = 1/2 
at singly occupied sites. Non-vanishing values of p^. and 

Xj_ indicate an EI state, which demonstrates the im- 
portance of the (transverse) 'spin-flip' processes for the 
formation and maintenance of f-c coherence. The slave- 
boson band shift \2Xz \ m Eq. (pT|) increases with in- 
creasing U (just as the Hartree shift). Obviously, the 
area of the EI phase is enlarged if one reduces the split- 
ting of the / and c band centers (cf. the red dot-dashed 
lines). We include the data for the metastable EI solution 
at E'^ = (as discussed above, in this case the charge- 
density- wave state will win), in order to show that (P 
and stay finite for all U. That means, for the orbital- 
degenerate EFKM (Ai^' = = 0), our 5'0(2)-invariant 
slave-boson scheme will not give the (artificial) transi- 
tion into an insulating Brinkmann-Rice-like correlated- 
insulator state'*^. This transition is a well-known short- 



FIG. 2: (Color online) [/-dependence of slave-boson fields and 
Langrange parameters for different Ef at T = {k — —0.8). 
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FIG. 3: (Color online) Band-renormalization factors at T = 
within S'0(2)-invariant (left-hand panel) and scalar (right- 
hand panel) slave-boson theory. Again, k — —0.8. 

coming of the scalar slave-boson approach to the Hub- 
bard model""^ and has been also observed applying the 
scalar slave-boson theory to the EFKM^. The effect be- 
comes even more apparent by comparing the variation of 
the slave-boson band-renormalization factors z^. 

Figure [3] shows that vanishes within the scalar slave- 
boson theory (right-hand panel) for iJ^ = at a critical 
interaction strength (C/br — 14.5), indicating the local- 
ization of charge carriers, whereas in our theory the band- 
width will be only slightly renormalized at this point (see 
left-hand panel) . Interestingly, the band renormalization 
is rather small in the EI phase as well (cf. the curves 
for = -1.2, -2.4). Here we find > 0.95, which 
explains the small deviation of the slave-boson order pa- 
rameter from its Hartree-Fock counterpart (see inset of 
Fig. p. 

Next we discuss the finite-temperature behavior. The 
variation of the EI order parameter and of the band- 
renormalization factor with T at fixed U is displayed 
in Fig. m Most important, in comparison with the 
Hartree-Fock data, the critical temperature for the 
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FIG. 4: (Color online) T-dependence of the EI order parame- 
ter A± and of the band renormalization z'^ at fixed Coulomb 
interaction (7 = 6 (left-hand panels) and U — 8 (right-hand 
panels). Red dotted lines show the corresponding Hartree- 
Fock data, where — 1. Band-structure parameters are the 
same as in Fig. [1] 




FIG. 5: T-dependence of the various slave-boson fields for 
U = 6, Ef = -2.4, and k = -0.8. 



El-semimetal/semiconductor phase transition is signifi- 
cantly reduced (see also Fig. [1}. Looking at z'^{T), this 
may be attributed to the more precise treatment of cor- 
relations and occupation number fluctuations. At Tc, the 
order parameter vanishes, and we observe a cusp in z^. 
Enhancing, above Tc, the temperature further, the band 
renormalization goes on, where now always decreases 
with increasing U. 

Figure [S] shows the temperature dependencies of the 
various slave-boson fields and Lagrange parameters for 
U — 6 (corresponding to the left-hand panel of Fig. 2]) . 

As expected, p± and are monotonously decreasing 

functions of T, with p±{Tc) = x'fiTc) = 0. The other 
fields exhibit a cusp structure at Tc- At higher tempera- 
tures the probability of finding double occupied sites and 
empty sites increases. At the same time, we find less 
singly occupied sites (cx (pq +pI)), which means that the 
increase of Pq is overcompensated by the reduction of p^ , 
indicating a more balanced occupation of / and c sites. 

Finally, we analyze the partial / and c electron density 
of states (DOS), p\{E) and Pi{E), defined via 



(81) 

(82) 
(83) 



where n^^^^^ are the corresponding particle densities. Fig- 
ure [S] gives p^^i{E) at the characteristic U-T points 
marked by LIV in the phase diagram of Fig. [1] Ob- 
viously, the high-temperature phase may be viewed as 
a metal/semimetal (panel I) or a small-gap semicon- 
ductor (panel II) in the weak-to-intermediate or strong 




FIG. 6: (Color online) Partial densities of states for /-band 
(red solid curves) and c-band (blue dashed curves) electrons 
at the points {U = 4.8, T = 0.45), (8.5, 0.45), (4.8, 0), (8.5, 
0) marked by I-IV in Fig. [T] Band-structure parameters are 
Ef ^ Ef = -2.4 and it = = -0.8. 



Coulomb-attraction regime. Accordingly, the EI phase at 
low temperatures shows different characteristics as well. 
As can be seen from panel III, a correlation-induced 'hy- 
bridization' gap opens in the DOS with n_ = 1 (n^ = 0) 
at T = 0, indicating EI long-range order. The pro- 
nounced c-f state mixing and strong enhancement of 
the DOS at the upper/lower valence/conducting band 
edges reminds a BCS-like pairing evolving from a (semi-) 
metallic state with a large Fermi surface above Tc- By 
contrast, the zero-temperature DOS shown in panel IV 
evolves from an already gapped high-temperature phase. 
Here, preformed pairs (excitons) may exist®'"^^, which un- 
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dergo a BEC transition at Tc- 

IV. SUMMARY 

In this work wc studied the extended Fahcov-Kimbah 
model with respect to the formation of an exciton con- 
densate, which is related to the problem of electronic fer- 
roelectricity. Motivated by the discrepancy concerning 
the existence of the excitonic insulator (EI) phase within 
the Hartree-Fock and scalar slave-boson approaches, we 
developed an S'0(2)-invariant slave-boson theory. The 
main result is that our improved slave-boson scheme is 
capable of describing the EI phase in a parameter region 
agreeing, at zero temperature, with Hartree-Fock (and, 
in 2D, constrained path Monte Carlo) results. This is in 
striking contrast to recent findings by the scalar slave- 
boson approach^l, which fails to detect the EI phase 
in the case of non-degenerate / and c orbitals. The 
agreement of the zero-temperature semimetal— >'EI and 
EI— )-band-insulator transition points with the Hartree- 



Fock and Monte Carlo values is ascribed to a rather 
weak band renormalization at T = 0. At finite tem- 
perature, band-renormalization effects due to electronic 
correlations and particle number fluctuations become im- 
portant, and, as a result, our slave-boson theory yields 
significantly lower transition temperatures than Hartree- 
Fock. From the analysis of the partial /, c, and quasipar- 
ticle densities of states, in the EI phase a crossover from 
a BCS-type condensate to a Bose-Einstein condensate of 
preformed excitons may be suggested. The results of our 
investigations may form the basis of forthcoming stud- 
ies, e.g., on the effects of fluctuations around the saddle 
point, allowing the calculation of pseudo-spin and charge 
susceptibilities for the EFKM on an equal footing. 
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